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Geodetic data at a plate boundary can reveal the pattern of 
subsurface displacements that accompany plate motion. We model 
these displacements as the sum of rigid block motion and the elastic 
effects of frictional interaction between blocks. We represent the 
frictional interactions by uniform dislocation on each of several 
rectangular fault patches. We then estimate the block veloci ties and 
fault parameters from geodetic data. Our Bayesian inversion 
procedure employs prior estimates based on geological and 
sei smologi cal data. We apply the method to the Transverse Ranges, 
using prior data from Bird and Rosenstock (1334) and geodetic data 
from the USGS trilateration networks. Our model consists of 11 
blocks and 26 rectangular fault patches. The block motion inferred 
from the geodetic data has the same order of magnitude as the 
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ABSTRACT 



Un i v e r s i 



geologic estimates, and for many faults the agreement is excellent. 
However , the geodetic data imply a displacement rate of about 20 
mm/yr across the San Andreas Fault, while the geologic estimates 
exceed 30 mm/yr. The prior model and the final estimates both imply 
about 10 mm/yr crustal shortening normal to the trend of the San 
Andreas Fault. Most of this shortening occurs on the Sierra 
Madr e-Cucamonga and the White Wolf Fault systems. Aseismic fault 
motion is a major contributor to plate motion, and the thickness of 
the frictional surface varies considerably from one fault to 
another. The geodetic data can help to identify faults that are 
suffering rapid stress accumulation; in the Transverse Ranges those 
faults are the San Andreas and the Santa Susana. 
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INTRODUCTION 


Since 1971 the U.3. Geological Survey (USGS) have carried out 
precise length measurements on baselines near the San Andreas Fault 
system in California (Savage et al . , 1381). These measurements tell 
us much about the details of plate motion, and the process of stress 
accumulation leading to earthquakes. 

In this paper we address the following questions: (1) How wide 
is the plate boundary, and can the plate motion be blamed on 
specific known faults? (2) Do the geodetic data agree with 
conclusions based on geologic observations and plate tectonic 
models? (3) Which faults are accumulating stress most rapidly? 

King and Savage (1334) analysed tr i later at i on data for the 
Transverse Ranges using a simple dislocation model. They showed that 
the strain rate is relatively low in this region, compared to other 
locations along the San Andreas Fault. This implies that the 
displacement rate on the San Andreas at depth is lower here than 
elsewhere, or that the San Andreas is locked to a great depth here, 
or possibly that displacement is taking place on faults outside of 
the tr i later at i on network. They found that there were significant 
spatial variations in strain rate, and that the data could be fit 
reasonably well by a model with only two faults: the San Andreas and 
the Gar lock. Their preferred model had 20 mrn/yr of right lateral 
slip on the San Andreas, and 3 mrn/yr left lateral slip on the 
Gar lock, with both faults locked to a depth of 15 km. 


Me are able to include many more faults than previous 
investigators because we use a new nonlinear inversion procedure 
incorporating prior estimates of all the parameters. The prior 
estimates are based on geologic and sei srnologi cal data. By 
including many mo re faults, we get a much more realistic model of 
the plate boundary region, and we can test the importance of many 
previously neglected faults. 

Our dislocation model is described more thoroughly in a 
separate paper (Matsu'ura et al . , 133S) reporting a similar analysis 
of the Hollister area of central California. The inversion method is 
described in detail in Jackson and Matsu'ura (1383). 



DISLOCATION MODEL 
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Assume t i o n s 


Me assume that geodetic displacements, and the block and fault 
motions that cause them, are constant in time over the period 
1371-1333. Thus we use rates of change of line length as our basic 
data, and velocities and displacement rates as our primary unknown 
parameters. By using geological estimates of block velocities in our 
prior model, we implicitely assume that displacement rates are 
constant over periods of many thousand years. The latter assumption 
is a working hypothesis, which we can test. If the final model 
agrees with the geological estimates, then the hypothesis of 
constant displacement rate cannot be rejected by our data. If in 
addition the parameters are well resolved by the geodetic data, then 
the final estimate is relatively independent of the prior model; in 
that case agreement between the prior and final estimates tends to 
confirm the constant rate hypothesis. 

Ne represent the crust by an elastic half -space divided into a 
finite number of blocks. The fault surfaces separating these blocks 
are divided somewhat arbitrarily into segments from 10 to 100 km 
long. In the absense of friction, the blocks would slide freely with 
no shear stress accumulation. However, friction on the upper part of 
the fault surface restricts motion and causes stress. Ne divide 
each fault segment into an upper "brittle zone" and a lower "ductile 
zone". The depth of the boundary between these zones 


i s cal 1 ed the 
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"locking depth;" the width of the brittle zone, measured in the 
fault plane, is called the "fault width." For vertical faults the 
locking depth and fault width are identical. In the ductile zone, 
the displacement rate across the fault is simply the velocity 
difference of the two blocks separated by the fault. In the brittle 
zone, the net displacement rate is the difference between the 
relative block velocity and the dislocation rate. This dislocation 
rate is introduced to represent the effects of frictional 

i n ter act i on j^and it is assumed to be constant over each dislocation 
patch. The dislocation rate is essentially a displacement deficit, 
likely to be repaid in the form of earthquakes or other episodic 
displacements at a later date. A schematic view of a single fault 
segment, and the notation we use to describe it, are shown in Fig. 
1 . 

If the entire surface separating two blocks were sliding 
freely, there would be no dislocation motion nor stress 

accumulation, and the surface displacement would be rigid block 
motion only. If the upper fault patch were completely locked, then 
the dislocation motion would equal the block motion, and stress 
would accumulate at a rate depending on the locking depth and the 
relative block velocity. 

We use the Jennings ej. al . (1375) fault map of California to 
determine block boundaries and the location of the upper corners of 
each fault segment. These are held fixed in our analy si s . . The width 
and dip angle of each fault patch, and the dislocation rate vector, 
are estimated from the tr i later at i on data. The block boundaries 
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that we assume are shown in Fig. 2, along with the actual traces of 

faults. 

% 

Over geologic time, the fitful motion in the brittle zone of 
any fault segment should average to the same rate as that of the 
lower ductile zone. Thus, the loading of the upper fault zone by 
frictional stress and its unloading by earthquakes are each 
temporary aberrations; the geological displacement rates should be 
compared to the steady "block motion" in our geodetic model. I«n the 
short term, friction at the upper part of the fault surface will 
cause some temporary distortion of the block: it is this distortion 

that we model with the dislocation. In order to distinguish the 
effects of block motion from the distortion near the block 
boundaries, we need to observe displacement both close to and far 
from the boundary . For many faults in the study area this condition 
is satisfied, and we can resolve well the dislocation motion that 
causes stress to build up on the fault. 

Calculating the rate of line length change caused by block 
motion is simple; the velocity of each monument is just the velocity 
of the block on which it rides. The rate of line length change 
between any two monuments is the projection of the velocity 
difference between the two sites onto the position vector from one 
site to the other. For the dislocations, we compute theoretical 
displacements and partial derivatives with respect to parameter 
values using the method of Matsu 'ura (1977), summarized briefly in 
Matsu 7 ur a ejt al . (1336). Me assume that the earth is a uniform 

elastic half-space, with the Lame' constants equal. 



P r i rn a r P a r am etevs 


The primary parameters consist of the eastward and northward 
velocity components for each block; and the dislocation rate (D), 
fault width (W) , dip angle (b), and slip angle (A) for each fault 
patch. In the Transverse Ranges we assume 11 blocks and 2b fault 
segments, so we have 22 block parameters and 104 fault parameters, 
for a total of 126. The dip angle is defined such that 90 degrees is 
vertical. The slip angle is- defined so that if the dislocation rate 
is positive, then 0 degrees represents left lateral strike slip 
motion, 90 degrees represents pure dip slip motion, and ISO A 
represents right lateral strike slip. If the dislocation rate is 
negative, then 0 and ISO degrees represent right and left lateral 
displacement deficits, respectively. 

Per i v e d P a r am e t e r s 

We also compute estimates of several derived parameters that 
are functions of the primary parameters. The "block slip" is the 
tangential component (parallel to the block boundary) of the 
relative velocity between two blocks. It depends on the block 
velocities and the orientation of the boundary. Positive block slip 
denotes right lateral motion. The "block convergence" is the 
component of relative motion normal to the boundary, measured such 
that convergence is positive. "Strike slip" and "dip slip" are the 
horizontal and updip components of the dislocation motion. A 
positive strike slip value indicates a right lateral displacement 
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deficit, and a positive dip slip value indicates a convergence 
deficit. "Creep rate" is the difference between the block slip and 

the strike slip. In other words, block slip is the tangential 
displacement rate at depth (below the dislocation patch), and creep 
rate is the shallow slip rate. 

Of course the earth is not a homogeneous half space, and fault 
displacement rates probably do not change di scon t i nuously at 
rectangular patch boundaries. The fault dislocation parameters 
represent averages over the fault zone, rather than specific values 
appropriate for any specific location. Similarly the "creep rate" 
represents the average shallow displacement inferred from geodetic 
baselines a few km long, and it might not be the same quantity that 
is observed with short baseline creeprneter s across some faults. 

For very long strike slip faults, the rate of stress 
accumulation can be calculated from the dislocation parameters; 
roughly, it is directly proportional to the dislocation rate and 
inversely proportional to the fault width. Assuming some fixed value 
for stress drop in a large earthquake, one can then calculate the 
recurrence timf and the characteristic displacement for such 
earthquakes. We give the relevant equations and calculate some 
estimated values of these parameters for faults in the Hollister 
Region in Matsu'ura e_t_ al . (1936). We do not report such 


calculations for 

the Transver 

se Ranges, 

because the 

required 

assumptions are ques 

ti enable in 

this region 

comp 1 i ca t ed 

by many 

fault intersections 

and d i p -si i p 

faults. 





INVERSION METHOD 


In our analysis the data are observed rates of change of line 
length, while the model parameters are the east and north components 
of each block velocity; and the dislocation rate, fault width, dip 
angle, and slip angle of each fault patch. The observed length rates 
are essentially linear functions of the block velocities and 
dislocation rates, but they are nonlinear functions of fault width, 
dip angle, and slip angle. Because of the large number of 
potentially active faults, the data may be insufficient to resolve 
all of ■ the relevant unknown parameters. Thus we are faced with a 
nonlinear, possibly underdetermined inverse problem. Howiever, we 
know a fair amount about the expected values of the parameters, 
independent of the geodetic data. For example, the block velocities 
should be the same order of magnitude as the geologically observed 
rates, and the dip and slip angle should agree with geological 
observations. Assuming that earthquakes occur in the brittle zone we 
can estimate the fault width at least approximately . It is then 
appropriate to use a nonlinear Bayesian procedure that makes use of 
the relevant prior information. Me use a method described in detail 
by Jackson and Matsu'ura (1335), and summarized only briefly here. 
Suppose that we have n observation equations and m unknown 
parameters, and let 

e = y - f ( x ) 

where y is an n -vector of observed data, x is an rn-vector of 
unknown parameters, f is a vector of possibly nonlinear functions 



giving the predicted values corresponding to the observations, and 
e is an h-vector of residuals. The equation above is referred to 
as the "observation equations." Suppose also that we have prior 
estimates x g for each of the parameters, and let 

d = x, - x 

be an m-vector of residuals to the prior estimates. If the errors in 
the observations and prior estimates can be described as random 
variables with mean equal to zero and covariance matrices E and 
D, respectively, then the minimum variance estimate of x 
minimizes 

i r -l r ^ 

T = e r E e + d D d 

Minimizing is equivalent to minimizing the sum of squared 

residuals to the combined set of equations 

e ' - F y - F f ( x) 
d' - Gx„ — Gx 

where F r E F = I and G T D G = I . The matrices F and G 
serve to standardize the observed and prior data, so that the 
equations above may be combined and solved by a standard nonlinear 
least squares estimation package . The prior estimates stabilize the 
inversion, but there is no need to treat them differently from the 


observation equations in computation. 
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Assuming that the functions f(x) are linear within a large 
enough neighborhood of the solution, the probable estimation errors 
are described by the asymptotic covariance matrix 

X = <£ t E‘ , A + of)' 1 

and the resolution matrix (Jackson and Matsu'ura, 1385) is 

R = X A r E H A 
= I 

R is an rn by rn matrix, whose diagonal elements give for each 
parameter the relative sensitivity of the final result to the 
observations. The complement (one minus the diagonal element) gives 
the relative sensitivity to the prior estimates. In the tables of 
parameter estimates below, the "resolution" for each parameter is 
the corresponding diagonal element of the resolution matrix, 
multiplied by 100 to convert to %. 

The reported error estimates for each parameter are the 
standard errors, or square roots of the relevant diagonal elements 
of the modified covariance matrix X', defined as 



where v L = T Z /n. '/is the "variance inflation factor," whose purpose 
is to adjust the uncertainties to match the observed sum of squared 
residals. We did not modify the covariances C and D, so that the 
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final error estimates for some parameters are larger than the prior 
uncertan.ties. » 
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T r i 1 a t e r a t i o n data 


We use line length data provided by Dr. James Savage and his 
group at U3GS in Menlo Park, California. For each line, we 
determined the rate of change of line length and its standard error 
by linear regression of length on time. The standard error is 
adjusted to be consistent with the length residuals, so that if the 
line length is quite linear with time, then the standard error will 
be small. In general, the standard error will decrease in proportion 
to the time span covered by the measurements. We assume that data 
errors are uncorrelated, so that the data covariance matrix C is a 
diagonal matrix whose diagonal elements are the squares of the 
corr espondi ng ' standard errors. We use data for 160 lines from the 
USGS Palmdale, Tehachapi, San Gabr i el , Los Padres, and San Fernando 
networks. The locations of the monuments are shown in Fig. 2. 'The 
length change rates and their standard errors are tabulated in Cheng 
(1335) and will be furnished on request. 

Pr i or est i mates 


We take prior estimates of block velocities from the model of 
Bird and Rosenstock (1334). Their model was adjusted to fit 
geologically observed displacement rates on major faults in southern 
California and the plate motion estimates of Minster and Jordan 
(1373). We make a few minor adjustments to their model: we modify 


some block boundaries slightly, and we combine their Santa Barbara 
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Channel, Chino, San Pedro, Santa Anna, and Called to blocks into a 
single block because they found negligible displacement between 
them. We assume rather generous prior uncertainties (20 rnrn/yr in 
each component) to allow for temporal changes from the geological 
average displacement rates. The prior estimates and their 
uncertainties appear in Table 1, along with our final estimates, 
which are discussed below. 

We assume prior estimates of the fault dislocation rate to 
offset the block displacements at the surface, except for a few 
faults such as "San Andreas F" where creep is observed. There, we 
choose a smaller dislocation rate, so that the difference between 
block slip and strike slip dislocation movement would equal the 
observed creep rate. We assume generous prior uncertainties for the 
dislocation rate as well; 20 rnrn/yr for most faults, and 10 rnrn/yr for 
a few less active faults. We assume that the fault width is 10 + 5 
km for all faults, based on the observation that most earthquakes , 
assumed to occur in the brittle zone, have depths in this range. 

A strike slip faults we assume to have a dip of 30 + 10 degrees, 

T»v 

except for the Santa Ynez and Pine Mtn faults which we give an 
uncertainty of 30 degrees. For thrust faults we assume a dip of '65 
+ 20 degrees. The Big Pine, San Cayetano, Sta Susana, San Gabriel, 
Sierra Madre, and Cucurnonga are assumed to dip down to the north, 
while the White Wolf and Pleito dip down to the south in our prior 
model. Prior estimates of slip angle are ISO degrees for strike slip 
faults, and 30 degrees for dip slip faults, with uncertainties 
ranging from 10 to 30 degrees depending on the geological complexity 
of the area. Table 2 shows the prior estimates and uncertainties of 
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all fault parameters, arid Table 3 shows the derived parameters 
corresponding to the prior estimates. 



RESULTS 




Goodness- o f f i t 

The block and fault model fits the trilateration data 
reasonably well, although the final rrns residual is about 1.4 times 
the estimated standard deviation of the data. This modest 
discrepancy might be caused by non tectonic deformation at some 
sites, by motion of neglected faults, by systematic errors in the 
data, or even by bad luck (that is, random coincidence). The 
presense of a few large residuals (two with absolute residual 
exceeding four standard deviations) suggests that there are some 
problems with a few of the data, especially from lines involving 
stations Tenhi, Pe2, and Tuj rrnl . Nevertheless the fit is good 
enough that there is no cause to doubt the overall data quality, or 
the estimated standard deviations of the data. 

Block v e 1 o c i ties and rn issing slip on t he San An dr eas 

Block velocities are shown in Table 1 and Fig. 3. The block 
velocities are resolved into parallel (block slip) and perpendicular 
(block convergence) components for each block boundary in Table 3. 
Fig. 4 shows the block slip, and Fig. 5 the block convergence. 
Referring to Table 1, we see that except for the Mojave block, which 
is fixed as a reference by a very strongly constrained prior 
estimate, the block velocities are resolved very well. They differ 
substantially from the prior estimates, in sharp contrast with our 
results for Hollister (Matsu'ura e_t_ al . . 1336), where the 
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estimated block velocities agree remarkably well with the geological 
estimates. In the Hollister region, our geodetic analysis showed 36 

mrri/yr of relative motion for blocks within the array; this compares 

with 34 rnrn/yr estimated from geological fault displacements, and 

leaves about 20 mrn/yr to occur outside the network, presumably to 

the west, if the total motion is to match the plate tectonic 

estimates of 56 rnm/yr (Minster and Jordan, 1373). In the Transverse 

Ranges, the geodetic data account for only about 20 mrn/yr of net 

motion across the geodetic array (that is, between the Malibu block 

(F) and the Mojave block (K). Following Bird and Rosenstock (1334) 

we assume that the Mojave block travels about 14 mrn/yr in the 

direction N11W with respect to the stable Worth American continent; 

then the Malibu block is moving 33 mrn/yr at M23W with respect to 

North America. Thus, the geodetic data fall about 25 rnrn/yr short of 

matching the predicted plate motion in the Transverse Ranges, 

compared to 20 rnrn/yr at Hollister . However at Hollister, the 

geodetic data match the geological observations, while in the 

Transverse ranges they don't. 

At Hollister, much of the shortfall between geodesy and plate 
tectonics is easily explained by additional displacements west of 
the • "network . In fact Hall (1331) and others estimate that 10 to 15 
rnm/yr right lateral displacement may have occurred on the San 
Sirneon-Hosgr i fault during the last 3 Myr . In the Transverse Ranges, 
a similar amount of displacement could be blamed on offshore faults, 
but there is still a direct conflict between the geodetic and 
geological estimates of motion on the San Andreas Fault system. Qur 
prior model, based on the geological compilation of Bird and 



Rosenstock, has the Malibu block moving almost 4? rnm/yr at N52W with 
II1S+) 

respect to the Moiave block. Our final model gives 20 rnm/yr at N40W. 


Possible explanations for the discrepancy between geodesy and 
geology in the Transverse Ranges are (a) the geodetic and geologic 
rates are not comparable because of temporal variations in 
displacement rate, (b) the geodetic estimates are in error, or (c) 
the geologic estimates err. 

We do not believe that time dependent block motion, hypothesis 
(a) above, provides the answer. While the near surface fault motions 
vary with time, the deeper block motions appear to be quite steady. 
This assertion is based on the close agreement between geodesy and 
plate tectonics in the Sal ton Trough (Savage e_t. al . . 1 373 ; Cheng, 
1985), and on the close agreement between geodesy and geology for 
Hollister (Hatsu'ura e_t al_. , 1936). But of course this evidence 
is circumstantial, and our assumption of constant block velocities 
over many thousands of years may fail. 

Errors in the geodetic model (hypothesis b) could result from 
erroneous data, or from a mistake in modeling. Data errors are very 
unlikely to be a serious problem in this analysis. While some 
systematic errors may cause annual or other short period variations 
(Jackson jet_ al . . 1933), such errors would have very little effect 
on the secular rates of line length change (Cheng, 1935). The 
t r i 1 at er a t i on data for the Transverse Ranges and for Hollister were 
collected by identical methods, and no discrepancy occurs at 



Hollister. A possible modeling error could result from our assigning 
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a prior fault width of 10 +5 km. The estimated fault width is 
poorly resolved (see discussion below), so that the final fault 
width estimate is strongly influenced by the prior. If the San 
Andreas Fault were locked to a much greater depth, then some of the 
resulting strain would occur outside of the network, and greater 
block motion would be required to match the observed displacements. 
We performed some calculations assuming much greater fault width on 
the San Andreas, and we found that the estimated block slip could 
match the geological observations if the fault is locked to 25 km or 
more (about the thickness of the crust). This seems much too deep to 
us, because earthquakes on this section of the San Andreas rarely 
exceed 15 km depth (Webb and Kan amor i, 1385) , and because the 
estimated fault width rarely exceeds 15 km in areas where it is well 
resolved (Cheng, 1385; Matsu'ura e_t_ al . . 1338). Nevertheless, we ' 
cannot completely reject the idea that the San Andreas is presently 
locked to the base of the crust within the Transverse Ranges. 

Errors in the geologic rates of Bird and Rosenstock^are not 

f 1 9 3 A- 1 

unthinkable. They constrained their model to agree with plate 
tectonics, possibly causing them to blame the San Andreas for motion 
that actually occurred on unknown faults. The Bird and Rosenstock 
model is supported by good independent data for fault "San Andreas 
A," north of the big bend (Clark et_ al_. , 1385; Dickinson e_t_ al . . 
1372; Si eh and Jahns, 1334), where 30 rnm/yr seems to be the minimum 
believable slip rate. However, on sections B-E , there is little 
data, and that published only in abstract form (Rust, 1332). Data 
for the southern section (F) suggest that 30 rnm/yr is an upper 
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limit. Humphries and Weldon (1934) suggest a much lower rate for the 
San Andreas, more consistent with our final estimates. The last 
word is yet to be written on this question, and it is conceivable 
that after further i nu es t i ga t i on the geological estimates will come 
into agreement with the geodetic values. 

In summary , we cannot rule out any of the listed hypotheses to 
explain the discrepancy between our estimated block velocities and 
the geological prior model. It could be- that the two are not 
comparable because block rates vary, or that we have badly 
underestimated the locking depth of the San Andreas, or that the 
geological estimates are off by 10-20 mrn/yr . 

Convergence across the Transverse Ranges 

Our results confirm the assertion of Bird and Rosenstock (1334) 
that substantial shortening occurs across the Transverse Ranges. 
Table 3 and Fig. 5 show about 10 mm/yr of block convergence across 
the Sierra Madr e-Cucamonga Fault system, wi th comparable values for 
the Northern San Andreas, Pleito, and White Wolf Faults. The 
estimated convergence is statistically significant at the 35 % 
confidence level for the San Gabriel W, Sierra Madre, Cucamonga, and 
Pleito Faults. For statistical significance, we use the approximate 
criterion that the absolute value of the parameter estimate should 
exceed twice the standard error. We did not calculate resolution 
estimates of derived parameters, but it is clear that the estimates- 
of block convergence come almost entirely from the trilateration 


data because the block convergence is derived from the well resolved 



block velocities, and because the final uncertainly is much smaller 
than the prior uncertainty. 


Humphries and Weldon (1334) objected to the convergence implied 
by Bird and Rosenstoc^because such convergence should cause massive 
crustal accumulation where only modest thickening is observed. It is 
not impossible that the convergence in our model results from 
systematic error causing an apparent secular decrease in line 
length. However, possible systematic errors have been quite 
exhaustively studied (Savage and Prescott, 1333; Jackson jat. al . . 
1333; Savage and Gu , 1335) and none has been identified that would 
cause a spurious secular dilatation. Another possible explanation, is 
that end effects from fault motion outside the array (on the central 
San Andreas, for example) would cause local contraction. However, 
this explanation is inadequate to explain the observed widespread 
convergence. We believe our estimates indicate true tectonic 
convergence, and that some explanation must be found for the missing 
crust. 


Dislocation rates 


The dislocation rates, listed in Table 2, 
are generally well resolved, but only for a few faults far e/ t h ey~) 
statistically significant at the 35% confidence level. The 

dislocation rates are projected into their strike slip and dip slip 
components in Table 3. The faults with significant strike slip 
dislocation rates are the Santa Susana E. , and the San Andreas A-E . 
None of the geodetic monuments in our study .were close to the 
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southernmost section of the San Andreas (section F) so the 
dislocation rate there is uncertain and poorly resolved. Dip slip 
fault displacement is statistically significant at the 35 % 

confidence level for only four faults: Santa Susan a E, San Gabriel 
H, Cucamonga, and Pleito (Table 3). The larger values ( for the 
latter three faults) are not geophy si cally reasonable, and we do not 
suggest that they be taken seriously. These erroneously high rates 
serve to adjust the motion of monuments very close to the faults 
from those predicted by the block motion. Because of the steep dip 
on the faults, the large dip slip displacements have a relatively 
small effect on the horizontal component of motion, which is what is 
actually measured. In specifying prior information, we did not 
adequately constraint' the dislocation motion for dip slip faults. 
The erroneous dip slip motion will have little effect on other 
parameters for the San Gabriel W and Cucamonga, because of their, 
very shallow depth. However, the Pleito is deeper, and the 
unreasonable estimate of the dip slip motion could cause a moderate 
error in the block motion for the Pleito Hills block. 

Fault w i d t h 

The fault width, listed in Table 2, is generally poorly 

resolved and seldom. di ffers from the prior estimate by more than 

twice the standard error . Only for the San Gabriel W and the Sierra 

Madre Fault does the resolution exceed SO %; both are dip slip 

faults with very small estimated widths. For two reasons we do not 
believe these fault widths are reliable, even if the asymptotic 
variance is small. First, we have found that the calculated length 



rates are strongly nonlinear functions of the fault width for 
shallow faults. Second, we have not provided adequate prior 
information for dip slip faults, so that the near -surf ace dip slip 
fault motion is poorly constrained in our model. This problem may 
also affect the depth estimates for shallow dip slip faults. 

In the Hollister area, Matsu' ur a e_t al ■ (193G) find 

resolvable depth variations on strike slip faults without the 
ambiguities present in the Transverse Ranges. Estimated fault widths 
varied from a few km to 14 km. As mentioned above, a fault width of 
25 km or more on the San Andreas would help to reconcile geodetic 
and geologic slip estimates, although the geodetic data cannot 
resolve this fault width. Nevertheless, the Hollister analysis 
shows that locking depth varies considerably from place to place, 
and the data for the Transverse Ranges add mild support. Local 
variations in earthquake depths (for example, Webb and Kan amor i, 
1935) add further support. 



CONCLUSIONS 


Estimated block velocities are well resolved. They show 
statistically significant strike slip at depth on the San Gabriel N, 
all segments of the San Andreas, the White Wolf E, and the Oarlock 
E. Cumulative strike slip motion is revealed in the relative block 
velocity between the Malibu and Mojave blocks, which amounts to only 
20 rnm/yr in the direction N40W. This estimate is approximately half 
of the geologically determined rate on the San Andreas. The apparent 
shortfall of geodetic slip could be explained by temporal 
variability of slip at depth, locking of the San Andreas to 25 km or 
more, or errors in the geological estimates. In any case there is 
geodetically observable slip on diverse faults, and the plate 
boundary region must be considerably wider than the. geodetic 
network, about 75 km in extent. The geodetic data show unambiguous 
crustal shortening in excess of 10 rnm/yr normal to the San Andreas 
system, mostly on the Sierra Madr e-Cucamonga and White Wolf Faults. 
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Table 3. Derived Fault Parameters 
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FIGURE CAPTIONS 


Fig. 1. Fault geometry and notation. The dislocation rate (D)., 
the fault width (W) , the dip angle («$"> , and the slip angle (X) are 
all treated as unknown parameters. The depth to the upper. fault edge 
(d) is fixed. 


Fig. 2. Map showing faults, both as mapped and as idealized 
for this study; blocks; and t r i 1 a t er a t i on monuments (triangles). 
Blocks are named in Table 1, and faults named in Table 2. The area 
shown is 210 km by 170 km. 

Fig. 3. Prior and final estimates of block motion. Final 
estimates indicate slower displacement across the San Andreas Fault 
than do the prior estimates. 

Fig, 4. Estimated block slip (that is, parallel slip at depth). 
Shaded arrows indicate right lateral motion, open arrows left 
lateral. 


Fig. 5. Estimated block convergence. Arrows pointing towards 
faults indicate convergence, arrows pointing away indicate 
di vergence . 










